Transfer and scattering of wave packets by a nonlinear trap 
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In the framework of a one-dimensional model with a tightly localized self-attractive nonlinearity, 
we study the formation and transfer (dragging) of a trapped mode by "nonlinear tweezers" , as well 
as the scattering of coherent linear wave packets on the stationary localized nonlinearity. The use of 
the nonlinear trap for the dragging allows one to pick up and transfer the relevant structures without 
grabbing surrounding "garbage" . A stability border for the dragged modes is identified by means of 
' of analytical estimates and systematic simulations. In the framework of the scattering problem, the 

shares of trapped, reflected, and transmitted wave fields are found. Quasi-Airy stationary modes 
^ , with a divergent norm, that may be dragged by the nonlinear trap moving at a constant acceleration, 

O ■ are briefly considered too. 
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C/^ i I. INTRODUCTION 

(~| . Solitons, in the form of robust localized pulses propagating in nonlinear dispersive media, have been studied exten- 
sively in various physical contexts Many works have been dealing with the soliton dynamics in media featuring 
uniform nonlinearities, modeled by wave equations with constant nonlinearity coefficients. More recently, much at- 
tention has been drawn to the study of solitons in structured media, which feature spatial modulations of the local 
nonlinearity strength (see recent review 0), including the case when stable bright solitons may be supported by 
self-defocusing nonlinearity whose strength grows towards the periphery W\. In solid-state physics, effective potentials 
induced by the nonlinearity modulation are often called pseudopotentials [A] . Recently developed experimental tech- 
niques make it possible to realize such structures in photonics and Bose-Einstein condensates (BECs). Particularly, in 
photonic media, the modulation profiles can be induced by nonuniform distributions of resonant dopants affecting the 
, local nonlinearity, which may be created by means of in-diffusion technology § • Another possibility is to use photonic 
crystals with voids infiltrated by liquids, which are matched to the host material in terms of the refractive index, but 
feature a mismatch in the Kerr coefficient j6i]. On the other hand, in the context of BECs, the local nonlinearity can 
be spatially modulated by means of the Feshbach resonance, induced by nonuniform magnetic or electric [8| fields, 
as well as by an appropriate optical field 

One-dimensional (ID) solitons were studied theoretically in the context of pseudopotential structures in BECs 
ITT| and optics The possible use of ID nonlinearity- modulation profiles was also elaborated in models of matter- 



13 , [IJ] (although the first experimentally 
15l|). Self- focusing pseudopotentials were 



wave lasers, designed so as to release periodic arrays of coherent pulses 
' realized prototypes of atomic-wave lasers were built in a different way 

theoretically studied in the two-dimensional (2D) geometry too, demonstrating that it is much more difficult (but 
possible) to stabilize 2D solitons in such settings, than using familiar linear potentials or lattices [iGl - fTsj . while the 
above-mentioned self-defocusing structures solves this problem easily [1]. 

The simplest example of the pseudopotential corresponds to the nonlinearity concentrated at a single point, which 
is accounted for by the (5- function |19|. In optics, it may represent a planar waveguide with a narrow nonlinear 
stripe embedded into it (stable spatial solitons supported by narrow stripes carrying the quadratic, rather than cubic, 
nonlinearity were found in Refs. 0, Hlj). In BECs, the localized nonlinearity may be imposed by a tightly focused 
laser beam inducing the Feshbach resonance. A relevant extension of such a nonlinear spot is the model of the double- 
well pseudopotential; in the simplest form, it is based on a symmetric set of two (5-functions multiplying the cubic 
term in the respective nonlinear Schrodinger (NLS) equation 22], or the quadratic nonlinearity in the model of the 
second-harmonic generation (2]| . In either case, full analytical solutions have been found for symmetric, asymmetric 
and antisymmetric localized modes trap ped by the set of the two ^-functions. In addition to that, full symmetric and 
asymmetric solutions can be found [23| for a similar discrete model, with two cubic-nonlinear sites embedded into a 
linear host chain, which was introduced (without the consideration of the spontaneous symmetry breaking) in Refs. 
[23 |. The symmetry-breaking problem was also recently analyzed in the two-dimensional setting, for a symmetric pair 
of nonlinear circles embedded into a linear host medium Il8i . 
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The model with the single nonlinear (5-function readily gives rise to a family of exact soliton solutions — see Eq. 
([5]) below — which are completely unstable, but may be stabilized if the point-like nonlinearity is combined with a 
periodic linear potential (in the same setting, the localized nonlinearity with the repulsive sign may support stable 
gap solitons) 'll]]. In the present work, we demonstrate that nonlinear self-trapped states may also be stabilized, 
without adding any linear potential, by the regularized localized nonlinearity, if the ideal (5-function is replaced by a 
Gaussian profile. 

In addition to supporting stationary trapped modes, the localized potentials may be used as tweezers, for the 



transfer of the trapped modes, which is a topic with many potential applications [25[. In particular, the use of the 
(5-functional linear trapping potentials makes it possible to find exact solutions for various transfer problems [26l . [27| . 
The dynamical process of extracting matter-wave pulses from a BEG reservoir by nonlinear tweezers, induced by an 
appropriate laser beam, was simulated in detail in Ref. [T^ . 

The main objective of the present work is to analyze possibilities of the controlled transfer of localized wave 
packets by a moving nonlinear trapping pseudopotential {nonlinear tweezers), represented by a narrow Gaussian 
profile multiplying the self-focusing cubic nonlinearity. An advantage of using the nonlinear tweezers is that they 
would not grab and drag small-amplitude (radiation) "garbage" , surrounding the target object in a realistic setting, 
and one may use the amplitude of the object to control its transfer (this resembles the well-known advantage of 
nonlinear optical amplifiers in comparison with their linear counterparts — see Ref. [1^ and references therein). The 
model is introduced in Section 2, where we also produce simple analytical results for the transfer problem, obtained 
by means of the adiabatic approximation, which is appropriate for a slowly moving trap. Basic numerical results are 
reported in Section 3. We also briefly consider a related problem of dragging wave patterns with a divergent norm 
by the nonlinear trap moving at a constant acceleration, which is relevant in connection to the recently considered 
transmission regimes for Airy beams, exact (29| or approximate (soj . including their nonlinear generalizations [3l|, 
and the transfer of linear trapped wave packets by a potential well moving at a constant acceleration [32] . In Section 
4, we report results for another natural problem related to the present setting, namely, the scattering of linear wave 
packets on the localized stationary attractive nonlinearity. Section 5 concludes the paper. 

II. THE MODEL AND ADIABATIC APPROXIMATION 
A. The formulation 

The model with the tightly localized nonlinear trap, which represents the moving tweezers, is based on the following 
normalized version of the NLS equation for wave function 'ip{x, t): 

fVt = -^^..-^~(x-e(i))iV'i'V', (1) 

where subscripts denote partial derivatives, while 5{x) and ^(t) define the shape of the nonlinear trap and its law 
of motion. In the context of BEGs, Eq. ([T]) is a normalized version of the Gross-Pitaevskii equation, with time t 
and coordinate x [note that in an early work [l9{. Eq. ([T]) was introduced as a model for tunneling of attractively 
interacting particles through a junction]. In the optical model, t is replaced by the propagation distance (2), while x 
is the transverse coordinate in the respective planar waveguide. 

In the first version of the model, the nonlinear trap was taken in the form of the ideal delta- function, 5 [x] = 5{x) 
[l9| . In the present work, numerical results are reported for its practically relevant Gaussian regularization, namely 

J(x) = (2^^i)-^xp(-xV(4e)), (2) 
with a sufficiently small width e. The law of motion for the moving trap will be taken as 

{0, at t < 0, 

\S 1 + sin {tt^-^^^ , at < t < T, (3) 
S at i > T, 

which implies that the trap is smoothly transferred from the initial position, ^ 0, starting at <: = 0, to the final 
position, ^ = S, at f = T. Both the initial and final velocities corresponding to this law of motion are zero, i.e., 
di/dt{t = 0) = di/dt{t = T) = 0. 
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B. Pulse solutions 

The model with the quiescent trap = 0), described by the ideal 5- function, allows one to construct solutions to 
Eq. ([T]) as a combination of two solutions of the linear Schrddinger equation in free space at a; < and a; > 0, coupled 
by the derivative-jump condition at a; = (while the wave function itself must be continuous at this point): 

= +0) - yj^ {x = -0) = -2 \yj{x - 0)1" yjix = 0). (4) 

This condition gives rise to a family of obvious solutions to Eq. ([T|) with d{x) — 5{x) and f = 0: 

where < is an arbitrary chemical potential. Such soliton-like solutions, with a discontinuous first derivative, are 
usually called peakons (see, e.g., Refs. |33|). Note that the norm of the peakon family does not depend on the chemical 
potential, namely, 

/ + 00 
|7/;(a;)pda; = 1. (6) 
-oo 

This degeneracy resembles the well-known feature of the Townes solitons in the two-dimensional NLS equation with the 
uniform nonlinearity [sS] . Accordingly, the application of the Vakhitov-Kolokolov ( VK) stability criterion, dN/dfi < 
[34. 35]. formally predicts the neutral stability of the peakon family; nevertheless, a numerical study demonstrates that 
the entire family is unstable [ll|, the instability being similar to that of the Townes solitons in two dimensions j34| . 
As mentioned above, the peakons can be stabilized by the addition of the linear periodic potential (with the nonlinear 
(5-function placed at an arbitrary position with respect to the potential (Tlj). In the next section, we demonstrate 
that stabilization may also be induced by the regularization of the i5-function, as per Eq. Generally speaking, 
it is easy to stabilize the family of Townes-like solitons, because the linearization of the underlying equation with 
respect to small perturbations around the soliton does not give rise to any unstable eigenvalue; in fact, the instability 
is nonlinear, i.e., it grows not exponentially, but rather as a power of time, and is represented by a respective pair of 
zero eigenvalues. Therefore, any modification of the equation which shifts the zero eigenvalues in the direction of real 
frequencies may stabilize the entire family. This will be the case in our considerations below (in Section 3) with the 
finite-width Gaussian ([2]) replacing the (5-function. 

Assuming that a stabilization mechanism is in action, but the shape of the soliton does not substantially deviate 
from Eq. ([5]), it is natural to expect that the trap, slowly moving according to given ^(i) [see, e.g., Eq. may 
drag the trapped soliton, which will then be described by the following modification of solution (O, in the adiabatic 
approximation: 

■0 = (— exp 

In fact, expression ([7]) is the usually defined Galilean boost of peakon ^ moving at instantaneous velocity d^/dt. 

C. Airy waves 

For the consideration of the nonlinear trap moving with constant acceleration W, we set ^{t) = {W/2)t'^. Then, 
Eq. ([1} can be rewritten in the accelerating reference frame moving along with the trap, i.e., in terms of the following 
variables: 

ij{x,t) = {W/A)^^^(l){z,T)exp[iWt{x-Wt'^/3)], (8) 
z = {2W)^/^x-Wty2), T={W^/2f'^t. (9) 
The accordingly transformed Eq. ([Ij reads: 

j0r+0^z-2<^ + 5(2-20) 101' = 0, (10) 

where zq is the location of the nonlinear trap with respect to the accelerating reference frame. Stationary solutions 
to Eq. (jlOp are 4'{z, r) = a(z), with real function a{z) obeying the following equation: 

-—^~za + 5{z-zo)a^=Q. (11) 



-iiit + - / — at + i—x 
^ 2 \dt dt 



exp 



(7) 
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Without the nonhnear term, Eq. (1111) is the classical Airy equation, whose relevant solutions at z < zq and z > zq 
are, respectively, given by: 



a-{z) = CiAi(z) + C2Bi(z), a+{z) = C3Ai{z). 



(12) 



Here, Ci and C2 are constants, and Ai(z), Bi(z) are the standard Airy functions which are defined by their asymptotic 
forms at z — >■ —00 and z — >■ +cxd 1361 : 



{Ai,Bi}(2 



1 



^i-z) 



— {sin, cos} 7 + ^ i-z) 



n3/2 



(13) 



{Ai,Bi} (z) 



20Fzi/' 



■exp 1 T-^z^^^ 



(14) 



Using the Wronskian of the Airy functions, which is equal to I/tt [33, the continuity condition for (j>{z) and jump 
relation 1^ at z — zq give coefficients Ci_2,3 [see Eq. (|12p ] in terms of the amplitude of the solution, A = a(z = zq), 
and Zq itself, that may be considered as two free parameters of the solution family: 



Ci = [(Ai(zo))" - 7rBi(zo)A2 
C2 = 7^Ai(zo)^^ 



A, 



C3 - AiAi{zo)r\ 

The known peculiarity of solutions based on the Airy functions is the divergence of the norm in the tail at z 
as the corresponding average density produced by Eqs. (IT^ and p^ . 



{a\z)) 



27rv^ ■ 



(15) 

—00, 

(16) 



decays too slowly. In time-dependent settings, the divergence of the norm implies a gradual decay of the initial pulse, 
which loses its power sucked into the tail growing toward z — )• —00 j30| . 



III. NUMERICAL RESULTS FOR THE TRANSFER AND DRAGGING PROBLEMS 



A. Stability conditions for the transfer 

First, we have constructed a family of stationary solutions to Eq. ([T]) with the quiescent (f — 0) regularized 
nonlinearity profile, taken as per Eq. ([2]), in the form of ip{x,t) — e~*^*u;(a;), where the real function w{x) satisfies 
the equation: 

l(fw ~. ,3 , 
^"^"'2^'^^^ (17) 

A typical example of the solution is displayed in Fig. [IJa), and the entire family of solutions is represented in 
Fig. [TJb) by the corresponding N{fi) curve, cf. the degenerate dependence given by Eq. © for the ideal 5-function. 
It is obvious that the entire family satisfies the VK criterion, dN/dji < 0, hence the solitons may be stable. This 
conjecture has been verified through the computation of the stability eigenvalues for perturbations around the steady- 
state solution of Eq. ((T7)) . and also by means of direct simulations of the perturbed evolution (not shown here). All 
the eigenvalues were found to be purely imaginary (they are defined so that imaginary values correspond to the real 
frequencies, i.e., neutral stability), therefore only the edge of the (continuous) spectrum j3^] lying on the imaginary 
axis is shown in Fig. [Ijc). This stabilization is a direct result of the regularization of the 5-function in Eq. p7)) . We 
have checked that, as smaller values of the regularization parameter, e, are taken, the corresponding stationary wave 
function w approaches the peakon profile of Eq. ([S]). It is worthy to note that, as the limit of the ideal (5- function is 
approached, we observe a pair of eigenvalues that bifurcate from the edge of the continuous spectrum and approach 
the origin, where they formally arrive in the limit of e — > 0, at which point the scaling invariance associated with 
Eq. ^ is restored, and the nonlinear instability emerges. 

We now turn to characteristic examples (cf. Fig. [5]) of the transfer of originally stable quiescent solitons by the 
nonlinear trap moving according to Eq. ([3]). As might be expected, the soliton survives the transfer, provided that it 
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is not dragged too hard, i.e., the transfer time, T, is not too short. We define the soliton as "surviving the transfer" 
if, at time t — T, its amphtude is no smaller than a sufficiently high fraction of the initial soliton amplitude — say, 
90% or 80%. For instance, in the case shown in Fig. [21 the soliton dragged by the trap of width e = 0.05 keeps 90% 
of the initial amplitude for T > 20. It can be checked that, in all cases of the "survival", the motion of the dragged 
soliton accurately follows the prediction of the adiabatic approximation [Eq. ([7])], as is clearly seen in the example 
displayed in Fig. [2jc), where the actual trajectory of the dragged mode closely follows the adiabatic trajectory ^(t). 
The surviving solitons essentially preserve their original shape of the smoothed peakon, cf. Figs. [IJa) andmd,e,f). 
The three latter figures demonstrate that the eventual shape of the transferred mode is nearly the same for three 
different values of the trap's width, e = 0.05, 0.1, and 0.025, i.e., the eventual results are robust, withstanding the 
variation of essential parameters of the scheme within wide limits. 

Results of the systematic simulations are summarized in Fig. [3J where the survival border is shown in the plane of 
the transfer parameters, S and T [see Eq. ([3])], for the fixed initial soliton and two different survival criteria, based on 
the demand of keeping the amplitude at the 90% or 80% level. It may again be concluded that the results are robust, 
as they only weakly depend on the particular choice of the criterion of the soliton's survival. The roughly parabolic 
shape of the border may be explained by the fact that (as discussed in the previous section) the gradual destruction 
of the dragged soliton is determined by the absolute value of the acceleration of the moving trap. In the case of the 
motion law ([3]), the acceleration is 

d^^ = -2r^"H"^^J' ^ ^ 

hence its average absolute value in the process of the dragging is (^\(P£_/dt'^\') ~ ttS/T^. Assuming that the soliton 
survives when the absolute value of the acceleration does not exceed a a certain critical level, the latter result predicts a 
parabolic shape of the survival border, T ^ To confirm the reasonable accuracy of the prediction, we have drawn 
a parabola with an empirically fitted coefficient, T — VTttS (the dashed line), to be compared to the numerically 
determined survival borders displayed in Fig. [31 



B. Quasi-Airy profiles dragged at a constant acceleration 

We have also briefly studied profiles of waves dragged at a constant acceleration, i.e., those built as per Eqs. 
and (|15l) . and their stability within the framework of Eq. (|T0)) . A typical example of the stationary wave attached 
to the accelerating nonlinear 5- function, i.e., a solution to Eq. (jTl]) . is displayed in Fig. [H Direct simulations of Eq. 
()10|) with this stationary wave used as the initial configuration were performed in order to examine its robustness. In 
particular, a{z) was constructed and an approximate local minimum of the quantity + C| [which determines the 
tail's amplitude, (a^(z)), as per Eq. (|16|) ] was identified for appropriate values of A and zq. The temporal evolution 
of the resulting configuration is shown in Fig. [4l 

A more detailed analysis of the quasi- Airy waves in the present nonlinear model will be a subject of a separate work. 
We note here that if we use small A (such as A = 0.00057 in the left column of Fig. [4]), the apparent discontinuity of 
the beam at Zq is barely visible; in this small-amplitude case, the solution appears to be more robust, and moves to 
the left, preserving its shape. If A is larger {A = 1.14 in the right panel), a significant fraction of the solution gets 




FIG. 1: (Color online) (a) An example of a numerically found stationary solution of Eq. (|17|l with fi = —1 and S{x) taken per 
Eq. ([2]). (b) The norm of the solutions, defined according to Eq. ([6]), vs. chemical potential /i. (c) The edge of the continuous 
spectrum lying on the imaginary axis [the full eigenvalue is defined as iXi, the continuous spectrum occupying the area above 
the border shown in (c)]. These results are obtained for a width of the moving trap of e = 0.05. 
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trapped by the delta- function at zqj while the larger remaining fraction still moves to the left, and a smaller fraction 
of the initial power disperses to the right of zq. 



IV. SCATTERING OF COHERENT WAVE PACKETS ON THE LOCALIZED NONLINEARITY 

The scattering problem in the framework of Eq. ([T]) with the ideal quiescent 5- function was first studied in Ref. [l^ , 
where the scattering of an incident plane wave was considered. It was found that the localized self- attractive non- 
linearity gave rise to an accordingly localized modulational instability (MI) of the incident wave, provided that its 
amplitude exceeded a certain threshold (minimum) value. Unlike the commonly known MI in the NLS equation 
with the uniform nonlinearity, the localized MI features a complex growth rate (i.e., it corresponds to the oscillatory 
instability) . 

Here, we aim to study the scattering from the regularized Gaussian nonlinear potential of incident wave packets in 
the form of coherent states. The latter represent the fundamental solutions of the linear Schrodinger equation in the 
free space, 

Here, B and b are real constants that determine, respectively, the amplitude and inverse width of the initial wave 
packet, A= {1 + bH'^y^''^ B and = VlF^^+~P are the same characteristics of the spreading packet; the norm of 
the packet is iV = -^/tt/ {2b)B^, while c is its velocity. In particular, in the case of small B^, the scattering problem 
may be considered in the Horn's approximation (see, e.g., Ref. [13]), by setting ip{x, t) = '!/'in(a;, t) + V'Born (a;, t) , where 



b {x - cty 

1 + ibt 



— t 
2 



(19) 




T=17 T=20 




XXX 

(d) (e) (f) 

FIG. 2: (Color online) (a,b,c): Density contour plots illustrating the dragging of the soliton, shown in Fig. [Ha), by the nonlinear 
trap with width e — 0.05, moving as per Eq. ((3|, for three different values of the dragging time, T. The black solid lines show 
the corresponding trajectory ^{t). Adopting the criterion that the surviving soliton must keep > 90% of its initial amplitude, 
we conclude that it survives the transfer for T > 20. Panel (d) displays the shape of the soliton in the case of T = 20, i.e., at 
the border of the survival, at time t — 30 (some time after the completion of the transfer). Panels (e) and (f) display the same 
as (d), but for the moving trap with widths e = 0.1 and e = 0.025, respectively. 
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% 5 10 15 20 



FIG. 3: (Color online) The top and bottom solid curves represent the boundary above which the initial soliton with chemical 
potential /i = — 1, shown in Fig.[lja), survives, keeping more than, respectively, 90% or 80% of the initial amplitude hy t — 30. 
The dashed curve depicts the analytical estimate (see text). 




(a) (b) 



FIG. 4: The temporal evolution of a{z) with zo — —3.86235, for the quasi-Airy waves dragged by the accelerating nonlinear 
trap with width e — 0.05. The left and right panels show, respectively, the stationary solution with A = 0.00057, and the 
unstable solution with A — 1.14. 

the first correction is to be found from the linear inhomogeneous equation, 

1 2 

i (V'BoriOf + 2 (V'Born)^^ = -^{x) |V'in| lAin- (20) 

Equation ()20|) can be solved by means of the Fourier transform. Here, we instead focus on direct simulations of the 
scattering problem in a wide parametric region in the space of (i?, 6), with the (5- function regularized according to 
Eq. ([2]). As illustrated by the examples displayed in Fig. [SJ the simulations demonstrate splitting of the incident 
wave packet into three parts — transmitted, trapped, and reflected ones. Naturally, the share of the trapped norm 
increases with B, i.e., with the enhancement of the nonlinearity; on the other hand, the dependence of the observed 
phenomenology on the inverse-width parameter b is relatively weak. It is also worthy to mention that, as seen in 
Fig.ISfd), the trapped pulse performs oscillations around the localized nonlinearity, while the transmitted beam splits 
into jets — see Figs. [5l[c,d). 

Results of the simulations of the scattering are summarized in Figs. |6] and [71 where we show the shares of the 
trapped, reflected, and transmitted norm (or power, in terms of the optics model), as functions of the amplitude (B) 
and inverse width (6) of the incident pulse for a fixed velocity, c — 20, and two different values of the scatterer's width, 
e = 0.05 and 0.1, respectively. The plots correspond to the moment of time t = 0.5, taken after the the incident pulse 
has passed the scatterer. Note that the sum of the three shares is, by construction, identically equal to 1. Analysis 
of the numerical data demonstrates that, in the limit of large amplitude (B), the transmission decreases while the 
reflection of the wave packet from the nonlinear barrier gets enhanced. 

Further, the comparison of Figs. [6] and [7] shows that the increase of the width leads to a gradual change of the 
scattering picture. It is quite natural that the wider scatterer (Fig. [7]) provides weaker reflection but stronger trapping 
of the incident waves. Generally, the solution of the scattering problem is more sensitive to the width e than the 
results reported above for the transfer problem, as in the latter case the wave packet had a chance to adjust itself to 
the particular shape of the nonlinear scatterer. 
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FIG. 5: Examples of the interaction of the incident coherent pulse (|19p with the localized nonlinear potential corresponding to 
the regularized profile for e = 0.05, at different values of the pulse's amplitude and velocity. 
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FIG. 6: Panels (a), (b), and (c) display, by means of gray-scale plots, the trapped, reflected, and transmitted shares of the 
norm of the incident wave packet, as functions of its amplitude (B) and inverse width (6), for fixed velocity c = 20 and width 
e = 0.05 of the nonlinear scatterer, at the moment of time t = 0.5. The trapped share is defined by the norm confined in 
interval —1 < a; < +1. 



V. CONCLUSION 



In this work, we have considered the ID Schrodinger model with the tightly localized nonlinearity, to study the 
following dynamical problems: the existence and stability of standing waves trapped by this nonlinear potential; the 
controllable transfer of the localized mode by the moving nonlinear trap (nonlinear tweezers); and the scattering 
of coherent wave packets on the stationary localized nonlinearity. By means of systematic simulations, the border 
of the effective stability has been identified for the transfer problem. Relative shares of the norm (power) of the 
trapped, reflected, and transmitted waves were found, as functions of the amplitude and width of the incident pulse, 
for the scattering problem. The outcome of the transfer is less sensitive to particular characteristics of the nonlinear 
potential trap (such as its width) than the results of the solution of the scattering problem. The quasi- Airy stationary 
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FIG. 7: The same as in Fig. [6] but but for a broader scatterer, with e — 0.10. 

modes, dragged by the nonlinear trap moving at a constant acceleration, were briefly considered too, and were 
found to be progressively less robust for beams with the increasing amplitude. Some results were obtained in the 
approximate analytical form, such as the adiabatic approximation for the slowly dragged localized mode, and a 
qualitative explanation of the quasi-parabolic shape of the stability border. 

This work may be naturally extended in other directions. It is especially interesting to consider the two-dimensional 
version of the transfer problem, using the corresponding localized nonlinear trap in the form of a circle (TgI . [Tsj and 
setting it in motion. A challenging issue in this case is the competition between the trapping of the two-dimensional 
wave packet, and its propensity to the intrinsic collapse. It may be also interesting to consider the scattering problem 
in one dimension in the case of a localized repulsive nonlinearity. Such studies will be reported in future publications. 
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